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We present new many-parameter families of strongly and symmetric hyperbolic formulations of 
Einstein's equations that include quite general algebraic and live gauge conditions for the lapse. 

The first system that we present has 30 variables and incorporates an algebraic relationship 
between the lapse and the determinant of the three metric that generalizes the densitized lapse 
prescription. The second system has 34 variables and uses a family of live gauges that generalizes 
the Bona-Masso slicing conditions. 

These systems have free parameters even after imposing hyperbolicity and are expected to be 
' useful in 3D numerical evolutions. We discuss under what conditions there are no superluminal 

f^*) ' characteristic speeds. 

o ' 

I. INTRODUCTION 

Q\ • The use of hyperbolic formulations of Einstein's equations in numerical simulations has several advantages (see jlj] 
for reviews). In particular this includes the ability of giving boundary conditions that are consistent with both the 
evolution equations and the constraints O. At the discrete level well posedness supplemented by consistency — in 
order to make sure one is solving Einstein's equations — implies that one can find discretizations that are numerically 
stable in the sense of the Lax theorem, a fact that is otherwise not necessarily true || . 

Well posedness guarantees that there is a bound in the growth of the solutions that is independent of the initial 
data, but it does allow this bound to grow with time, and even fast. If this happens, numerical solutions can grow 
■ fast as well and, since Einstein's equations are non-linear, this can make the code crash in a finite time. This is what 
usually happens in 3D black hole simulations. This growth might appear due to several factors, including fast growing 
s ! . gauge or constraint violating modes. 

Initially, hyperbolic formulations of Einstein's equations with associated well posed initial-value problems relied 
on the use of harmonic coordinates Q, that is, V Q V Q a; M = 0. Although these coordinates have been used recently 
(see H and references therein, and also ||), they might be too restrictive for numerical applications, since one might 
^ | want to use another gauge which is better suited to the dynamics of a given numerical evolution. Initial efforts 
towards relaxing the harmonic condition introduced hyperbolic formulations that used time harmonic coordinates, 
i.e. V a V Q t = (sec & |tJ and references therein). This condition can be written as 
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■ where K is the trace of the extrinsic curvature, (3 l the shift vector, and N the lapse. A generalization of the time 



d t N-f3 l d l N = -N 2 f(N)K (1) 

with / > was introduced by Bona and Masso (BM) || (^o|, being able to incorporate this condition in a strongly 
hyperbolic (SH) evolution system. Here by SH we mean that the principal part has a complete set of eigenvectors with 
real eigenvalues. If the principal part can be diagonalized with a transformation that is uniformly bounded and smooth 
in all of its arguments, then the initial value problem can be shown to be well posed (this will be discussed later). 
The BM formulation represents the first effort in moving away from the time harmonic slice while achieving strong 
hyperbolicity, incorporating other gauge conditions often used in numerical relativity. For example, (Q) includes the 
"1 + log" condition (/ = a/N, with a some constant). In fact, given a hyperbolic formulation of Einstein's equations 
that introduces the lapse as a dynamical variable through condition (|l|), one can add to the right hand side (RHS) of 
this equation any function S of spacetime, 

d t N - P&N = -N 2 f(N)K + S(a?) , (2) 
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without changing the level of hyperbolicity of the whole system. From this observation it is clear that any spacetime 
with given lapse and shift satisfies equation (||), provided S is chosen appropriately (namely, S — d t N — f3 l diN + 
N 2 f(N)K). ' 

A closely related slicing condition, which has also been used in several hyperbolic formulations (see [7], [10| and 
references therein), is obtained by densitizing the lapse. That is, one writes 

N = e V , (3) 

where Q = Q{x^) is an arbitrary but a priori specified function of spacetime, g is the determinant of the three-metric 
and a any constant (strong hyperbolicity implies a > 0). As before, given any spacetime with lapse and shift there 
is always a Q that satisfies eq. (||) . Restricting to the case of zero shift and taking a time derivative of both sides of 
Eq.(||) one gets 

d t N = -2aN 2 K + Nd t Q , 

which shows that in this case condition (^) is time harmonic if a = 1/2 and dtQ = 0. 

Most hyperbolic formulations use either a densitized lapse or some of the BM conditions. The Kidder-Scheel- 
Teukolsky (KST) formulation |l0| ], for example, uses densitized lapse. This formulation is a many-parameter general- 
ization of previous systems ]TT| that is SH or symmetric hyperbolic if these free parameters satisfy certain inequalities. 
KST have shown that one can make use of the freedom in these parameters to extend the lifetime of 3D numerical 
evolutions of black holes. 

As already discussed, one can always evolve a known spacetime in a given slicing by using the appropriate den- 
sitization Q. Sometimes one can even use the same function Q to evolve a high, non-linear distortion of the same 
spacetime |l2[. But it is not clear what densitization to choose if one aims to evolve a spacetime that describes, say, 
the collision of two black holes. This is a rather general problem of gauge prescription. For example, even though 
there are some features of the BM slicings that are already understood (like its singularity avoiding properties), it is 
still far from clear what the "most appropriate choice" for the function f{N) is (see section [v|). However, it seems 
one could benefit from more general slicing conditions in the same way the free parameters of the KST system turn 
out to be useful even though one does not completely understand why. 

The aim of this paper is to present families of SH and symmetric formulations of Einstein's equations that combine 
both the freedom of the KST system with a family of lapse choices that includes those described above. 

The first family that we present is a straightforward generalization of the KST one that has, instead of a densitized 
lapse, an algebraic relation of the form 

N = N(jg,x*). (4) 

That is, the lapse still depends on the spacetime coordinates and the determinant of the three-metric in an a priori 
prescribed way, but now the power law dependency on g is relaxed. The motivation for introducing this generalization 
is that it is closely related to the BM conditions, in the same way densitized lapse is related to a time harmonic slicing. 
More specifically, condition (EJ) is equivalent to a particular case of condition (||) presented below provided the function 
F in Eq.(||) satisfies 8k F > or, equivalently, d g N > 0. As we will see later, these two conditions are necessary 
for strong hyperbolicity. Therefore, in those cases there will be a correspondence between the live gauges and their 
algebraic counterpart (though the resulting evolution equations will not necessarily be equivalent off the constraint 
surface). This family has, as the KST one, 30 variables and is presented in Section 0. There we also perform a 
characteristic analysis in order to give sufficient conditions for strong hyperbolicity to hold and in order to analyze 
under what conditions there are no superluminal speeds. 



In section III we present a many-parameter family of SH formulations that includes the following live gauge condi- 
tions 

d t N - (3 l d t N = -NF(N,K,x^) , (5) 

where F is any function of its arguments that satisfies djcF > 0. These are generalizations of the BM ones that 
relax the linear dependency on K. Since now there are evolution equations for the lapse and its spatial derivatives, 
our system has 34 variables. This system differs from the BM formulation in various aspects. For example, it uses 
different variables (in fact, the BM formulation introduces three extra variables in addition to the ones we consider 
here). 



In section IV we will explicitly show that some of the formulations presented in this paper are symmetric hyperbolic; 
well posedness of the initial value problem follows immediately for these cases. These symmetric hyperbolic systems 
are obtained by restricting some of the free parameters of our formulations, but not the gauge choice. That is, we are 
able to get symmetric hyperbolic systems with lapses given by Eqs. (Eh or (H). 
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II. STRONGLY HYPERBOLIC FORMULATIONS WITH AN ALGEBRAIC GAUGE 

The system presented in this section consists of the 30 variables {gij, Kij, dkij}, where g^ is the three-metric, 
the extrinsic curvature, and where the extra variables dkij are equal to the first order spatial derivatives d^g^ of the 
three-metric provided the constraints are satisfied. 

The evolution equations in vacuum are obtained by adding constraints to the RHS of the evolution equations 
obtained from setting to zero the four-dimensional Ricci tensor. Following the notation of KST (except that here we 
define d = {d t - £p)/N), 

do9ij = ~2Kij i (6) 

doKij = R l3 - iViV.-iV - 2K ia K a j + KK l3 +19lJ C + C g ab C a{ij)b , (7) 



-9.! 
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d d ki j = -2d k Kij - + r]g k ^C^ +X9ijC k , (8) 



where {7, £, 77, x} are free parameters, C = (R — K ab K ab + K 2 )/2 is the Hamiltonian constraint, Cj = V a K a i — ViK 
the momentum one, and C k ij = dkij — d k gij, Cikij — dud k ]ij are constraints that arise due to the introduction of the 
extra variables. The Ricci tensor R\j belonging to the three-metric is written as 

Rij = 2# ab {-dadbtj + d a d(ij)b + 9(id| Q (,|j) - d(idj} ab ) + -d i ab d jab + -{d k - 2b k )T k lj - V k ^T l ik , 

where bj = d k ijg kl , dk = dkijg 1 ^ and 

1 



Similarly, the momentum constraint gives 



d = g ab {d a K u - diK ab ) + \(d k - 2b k )K kl + ^ ab K ab . 

The shift is assumed to be a prescribed function of spacctimc. In contrast to KST, where the lapse is defined by 
Eq. (^), here we will consider the lapse to be an arbitrary function of the coordinates and the determinant of the 
three metric, N = N{g,x^) (as described below, strong hyperbolicity requires d g N > 0). In this case, we have 

d k N , 8N 

Cp.f f U k 



N ejJ K N8x k ' 

1 1 fd 2 N d 2 N d 2 N d 2 N 

7 ViVjN = <Teff(9 ab d (i d j)ab -d i ab dj ab + d i dj) + —(-^g 2 didj+g—— 1 d j +g- 1 



where the "effective" a is defined by 



ON 



{(Jeff coincides with a in Eq.Q if the lapse is densitized). In order to analyze hyperbolicity one has to look at the 
principal part of the system, that is, the terms that have spatial derivatives of the main variables. In this case the 
principal part is 

d Q gij = l.o. , (9) 

doKij = ^g ab {-d a d blj + (1 + Qd a d {i j) b + (1 - C)%^| a 6|j) - (1 + 2a V/) d {l d j)ab + 7 g l3 g kl d a {d klb ~ d bk ij) + l.o. ,(10) 

dodkij = -2d k K lj +i 1 gk {i g ab {d\ a \K j}b -d j) K ab )+xg lJ g ab {daKk b -d k K ab ) + l.o., (11) 

where l.o. stands for "lower order terms". The characteristic speeds in the direction n l are given by (3 % rii, ±7V + /3*rij, 
±Ny/\i + j3 % rii (see next subsection for a derivation and details), where 

Ai = 2a eff , 

A 2 = l + x-i(l + Ch + 7(2-?7 + 2x) , 

A3 = lx+l^-0v-\{l + 2o eff ){r, + 3x). 
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The system is SH if 

Aj >0, for j = 1,2, 3, 

A 3 = ^(3Ai + l) if Ai = A 2 . 

The system has no superluminal speeds provided that < A^ < 1 . One way of achieving this is by asking all the A.; 's to 
be or 1. Since Ai > for strong hyperbolicity, one needs Aj = 1. In particular, the condition Ai = 1 is possible only if 
the lapse is densitized with a constant <J e ff equal to 1/2. In that case one has two families of formulations with speeds 
along the light cone or normal to the hypersurfaces, one of them bi-parametric and the other one mono-parametric, 
see |10[. If (J e ff is not constant, one can ask all the speeds but Ai to be one. If one requires the parameters 7, 77, \ 
to be constant this leads to two mono-parametric families of SH systems: 

„ , 5 6 2 

a eff > °> C = --,7 arbitrary ,r) = -,X = —= 

or 

(Jeff > 0, C= ~ 5X n + 8 ,7 = -% r ) = -3x,X 7^ arbitrary 
9x 3 

As shown below, the characteristic modes (and whether they are superluminal or not) associated with Ai are inde- 
pendent of the formulation, they depend only on the choice of the slicing condition. 




A. Characteristic analysis 

Here, we discuss under which conditions the system (|^-[[l]) is SH, and what the characteristic speeds are. In order 
to do so, we choose a fixed direction n % with gijrtni = 1 and study the eigenvalues of the principal part of d^-pl) in 
the direction of n k : 

\ / 9ij \ 

A , (12) 

B J \ d kij J 

where the matrices A and B are obtained from the principal parts of ([l0|) and (|ll|), respectively, by replacing dk by 
rife . The characteristic speeds of the system are given by 

N^l + ftm , 

where /i are the eigenvalues of problem (|l2|). The system is SH if all the eigenvalues are real and the corresponding 
eigenvectors are complete. 

These eigenvalues are either zero or can be obtained by considering the eigenvalue problem 

^K,, = ABKij . (13) 

Explicitly, we have 

^ 2 K Z j = K i: j + Ari{in s Kj) s + BmrijK + Cgij (n r n s K rs — K) , 
where the coefficients A, B and C are 

A = -2 + X ~ 7(C - 1)'7 - (veff + 1)(V + 3x) , 



B = -X + -(C - 1>? + ke// + -) (2 + „ + 3x), 

-2C = X-^(C+l)»? + 7(2-r? + 2x). 

Next, we complete n % to a complex orthonormal basis n l , m', fh % such that n'"n J gij — 1, m % fh? g.ij — 1, and all other 
scalar products are zero. We then decompose K^ according to 

Kij = ariiiij + bm^ifhj) + [cn^mjj + dm^im^ + c.c] , 
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where a, b are real and c, d are complex. In this basis the linear operator AB takes the simple form 



AB 



The eigenvalues of this matrix are 



/ 1 + A + B B-C 
1 - 2C 

01 + 4 

\ 







Ai 


= 1 


+ A + 


A 2 


= 1 


-2C-- 


A 3 




A 


= 1 


+ 2 = 


A 4 


= 1, 





B = 2a, 



eff, 








(14) 



1 + X- o(l + 0»7 + 7(2-»7 + 2x), 



2*- 



|(1 - C)»? - \(l + 2a)(rj + 3 X ), 



and we have 



B — 



C = \ 1 + -(\ 2 + l)-2\ 3 . 



We demand that the matrix AB be diagonalizable and have only positive real eigenvalues. As we show now, this 
is a sufficient condition for the original system (|^-[ll[) to be SH. 

First note that since AB is positive definite, B is injective and A is surjective. Now let 
six real eigenvectors of AB which have the eigenvalues o;W, ...,w^'. Then, the 12 vectors 



K^' denote the 



V± s 





(Bif (s) )«j 



(15) 



1, 2, 6. These vectors are linearly independent since 



are eigenvectors of (|T^) with eigenvalues fi± s - 

c<j( s ' ^ and since B is injective. The remaining eigenvectors have \x = 0. Since A is surjective, we have 

dim ker A = 18 - dimQ?A = 12, 

thus there are 18 zero eigenvectors, 12 with non-trivial dki^s which lie in the kernel of A, and 6 with non-trivial (u, -'s. 
Therefore, we have a set of 30 independent eigenvectors with real eigenvalues. In order to show that the system (pHnf) 
yields a well posed initial- value problem, one has to look at the matrix S(n l ) whose columns are the 30 eigenvectors 
of the principal part of the system and show that S(n l ) and its inverse are uniformly bounded and that they depend 
smoothly on ri 1 and the metric coefficients gij and g 13 . We do not show this. However, if we linearize the equations 
around flat spacetime in Cartesian coordinates, the principal part depends only on ri 1 and the flat metric Sij. Using 
an isotropy argument it is not difficult to show that in this case the matrix S(n l ) can be obtained from S(n l ) in a 
fixed direction n l Q by a rotation which maps n l Q to n l . It follows from this that S(n l ) can be chosen such that its 



norm and the norm of its inverse are uniformly bounded for all ri 1 with Si 



guarantee well posedness 1 13 



1. In this case, this is sufficient to 



The matrix AB is diagonalizable if and only if Ai, A 2 , A 3 > and A 3 = (3Ai + l)/4 whenever Ai = A 2 . 



B. Characterization of the eigenmodes 

Here we give an interpretation to the eigenmodes of (|l3|) when the SH field equations (HJ7|^) are linearized around 
Minkowski spacetime. One can see that in such case the characteristic modes are solutions to Einstein's equations. 
According to (|lj) these modes are given by the following vectors corresponding to the values of (a, 6, c, d): 

A = l: (0,0, 0,1), (0,0,0,*), 

A = A 3 : (0,0,1,0),(0,0,*,0), 

A = A 2 : (0,1,0,0) if A 2 = A x , (C - B, 2a eff - 1 + 2C, 0,0) otherwise , 

A = Ax: (1,0,0,0). (16) 
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The first two modes, which propagate along the light cone (ie. which have speeds ±JV + /? l rij) are physical modes: 
By making a Fourier transform of the linearized Hamiltonian and momentum constraints, it is easy to check that the 
constraints are satisfied by the corresponding 4 eigenmodes (|l5|). Next, we have six constraint violating modes which 
have characteristic speeds ±Ns/X3 + (3 l rii, ±N^/X2 + {Fnc Indeed, the Fourier transformed linearized momentum 
constraint yields 

•n? Kij — riiK = —brii H — crrii H — crhi , 



and only modes with b = c = satisfy these constraints. The modes with characteristic speeds ±N \j2a e f / + (3 % rii are 
gauge modes: With respect to an infinitesimal coordinate transformation of the form Sx^ i— > Sx^ + fS^ , the extrinsic 
curvature transforms according to 

i-> K^ + ii.ii ,f , 

therefore 

an a + / , 

and the eigenmode (|l^) with eigenvalue Ai = 2cr e // can be gauged away. 



III. STRONGLY HYPERBOLIC FORMULATIONS WITH A LIVE GAUGE 



Here we construct SH formulations with live gauges. In addition to the variables of the previous section, we promote 
the lapse N and three extra quantities Ai which are equal to (diN)/N if the constraints are satisfied, to independent 
variables. Our variables are, therefore, {gij, K^, dkij , N, Ai}. As in the previous section, the shift is assumed to be 
an arbitrary but apriori prescribed function of spacetime. 

As evolution equation for the lapse we consider 

d N = —F(N, K, x 11 ), 

where F(N, K, x^) is an arbitrary function of its arguments. An evolution equation for A4 is obtained from this by 
taking a spatial derivative: 

9F . 1 dF n T ^ 1 dF 

doA t = — - — Ai ^—d l K — +£C l , 

dN NdK N dx l 

where we have also added the momentum constraint with a free parameter £. 
The evolution equations for g^, and dkij are the same as in (|||^) where now 

^ViV.-iV = d {l A 3) - T%A k + AiAj , 

and where we replace (dkN) /N by A^ in the evolution equation for dkij ■ 
The principal part of the system is 



d 9ij 
d K i:j 



l.o. 



g ab (-d a dbij + (1 + Qd a d {lj)b + (1 - (,)d(id\ ab \j) - d {l d o)ab + 7 gijg kl d a (d k i b - d bk i)) - <%A,) + l.o. , 



dodkij = -2dkKij +r)g k ( l g a {d\ a \Kj) b - dj)K ab ) + xgijg a (d a K k b - dkK ab ) + l.o 
d N = l.o., 



d Ai 



1 dF 

NdK 



8iK + Zg ab (d a K bi - diK ab ) + l.o. 



In order to get the conditions under which the system is SH, and in order to get the characteristic speeds, we can 
use the techniques of the previous section: The principal part has the form 




/ 



A 

Vo B 



/ N 

9ij 



K, 



1.1 
u 
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where u = (dkij,^i) ■ Explicitly writing the matrix AB, we obtain the same result as in ( |14| ) except that now 

A = -2 - |x- |(3C - - ^ 

B = ^(3C-l)»7 + £ + l + 2<re// 

-2C = X ~^(C + l)»? + 7(2 -»? + 2x), 

where er e // = (djcF)/(2N). Therefore, the characteristic speeds of the system are (3 l rn, ±N + f3 l rii, ±Ny/Xi + (3 l n i: 
with 

Ai = 2o- e// , 

A 2 = l + x-i(l + C)^ + 7(2-r? + 2 X ), 

a 3 = -\x-\(H-i)v-li- 

In particular, if one chooses £ = cr e // (tj + 3x) , these speeds are exactly those of the system of the previous section and 
strong hyperbolicity holds under the same conditions. The conditions under which there are no superluminal speeds 
are also the same, and the mode associated with Ai is related to the choice of gauge as well. More generally, one can, 
as in the previous section, ask all the Aj's but Ai to be 1. This leads to two SH many-parameter families, one with 
three free parameters (7,^,77): 



(Jeff > 0, 



X 



2 ' 

(1 + 077-27(2-7?) 
2(1 + 2 7 ) 



£ = ~X-^(3C-l)t7-2, 
and another one with two free parameters (C,x) : 

cr e// >0, 7=~2> C?7 = -2, £ = -\x+\ r l-\- 

IV. SOME SYMMETRIC HYPERBOLIC SUBFAMILIES 

Here wc show that some of our formulations are symmetric hyperbolic. We do not intend to give the most general 
conditions under which this holds but, instead, show that there arc at least some subfamilies that are symmetric 
hyperbolic. We start discussing the family of live gauge conditions and then briefly discuss the algebraic case. 

We show that the principal part of the system can be brought into symmetric form by using a transformation which 
does not depend on 71*. In order to find such a transformation, it is convenient to first transform the variables Kij 
and dkij into their trace and trace-less parts: 

= Pij + - gtjK, 

3 11 

dkij — "Zeuj + g 9k{iFj) — g ffijlfc + - g t jdk , 

where Tk — bk — dfc/3 and Pij and ekij are trace-less in all their indices. In terms of the new variables K, Ay, ekij, 
T/f, dk, the principal part is 

I^K = ( 1 + | 7 ) ( r„ - %dn) - A n , 



2 7V 3 

(iPij = -e-nij + (1 + e (tj) 



1 



— (5 - 9C)n(iTj) - - n {i dj) - n (i A 
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V - 2 )Pk 



10 



■qn k K, 



fid k = (rj + 3x)P kn - g (3 + f] + 3x)n k K 7 



liA k = £P, 



kn 



2<J eff + ) n k K, 



3 1 

-n k Pij + — g k (iPj) n — — gijP k n , 



where TF indicates the trace-free part, and where T n = n k T k , e n ij — n k ekij etc. First, we see that if we set £ = — 1, 
the term involving e k ij in the equation for P^ is the symmetric counterpart of the term involving Pij in the equation 
for ekij- Now, we can try to find three independent linear combinations of the variables T k , d k and A k and rescale K 
such that the principal part becomes symmetric. For the choice £ — — 1, 7 = —2/3, 77 = 6m/5, \ — — 2m/5, we find 
that the transformation 



/ 2m 

K^c x K, r fe H^c 2 (r fe —d k 

yields a symmetric system, provided that 



d k 



C3 



a eff + g£ ) d k - A k 



m > 1, cr e // > p = 



7m 
15" 



1 

6' 



and if we choose 



A k 1 ^ c 4 ^4 fe 



(1 + 3cla eff ) + J(1- H°eff) 2 + 12c|p 



2c| ( ^e// + 3^ 



4 



20(m- 1) ' 
1 + 



Note that the definition of £ and the requirement <r e // > p guarantee that — 1 < c|£ < and 3cr e // + £ > 0. The 
characteristic speeds of the system are P l rn, ±N + (3 l n>i, ±N\f\i + (3 l rii, with 



Ai = 2<7 e/ /, A 2 = 2p, 



7m I 

h To ~ 2 e 



Setting, for example, m = 15/14 and c| = 2 yields A 2 = 2/3 and A 3 = 3/4 — ^/2 < 1. Therefore, as long as <r e // > 1/3, 
the system is symmetric hyperbolic and has no supcrluminal constraint violating modes. 

One can similarly show that there are symmetric hyperbolic subfamilies in the algebraic case. For example, the 
choice 



c = -i, 



7 



V= y(l + 6cr e //), 



X 



y(l + 6cr e// ), 



yields a symmctrizable system provided that & e ff > 3/10. In this case the A^'s are Ai = A 2 = 2cr e //, A3 = (l+6cr e //)/4. 

Finally, as in KST, one could perform a many-parameter change of variables in any of our formulations without 
affecting the spectrum of the principal part and therefore without changing the level of hyperbolicity. 



V. DISCUSSION 



The freedom in the KST formulations has proven to be very useful in improving the stability of 3D single black hole 
numerical evolutions. There are studies underway to understand the reasons behind this, but there is still not a clear 
picture. If the lifetime of present evolutions using this system carry over to dynamical situations, interesting parts of 
a binary black hole collision could be described. However, one possible obstacle is the lack of flexibility of the family 
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of lapses considered in the KST formulation, namely the densitized lapse prescription. On the other hand, the main 
idea behind the BM formulation is to introduce in SH formulations several of the dynamical slicings that are used in 
numerical evolutions. The purpose of this paper has been to combine the spirits of the KST and BM formulations, 
and the result is expected to be useful in 3D evolutions, especially in dynamical ones. 

We have considered an algebraic and a live family of slicing conditions that include most of the conditions used 
in numerical relativity, such as densitized lapse, time harmonic slicings, the "1 + log" case, and the BM conditions. 
Furthermore, we have shown that one can obtain symmetric hyperbolic systems with these choices of gauges and, 
therefore, have a well posed initial-value problem. Up to our knowledge, this is the first time this has been achieved. 
The BM formulation, for example, is shown in |^| to be SH in the sense used in this paper but, as already mentioned, 
this does not automatically imply well posedness. 

Initially the approach to black hole evolutions (see 1 14 for a review in numerical relativity) was through the use 
of singularity avoiding (SA) slicings, such as maximal slicing. This gauge condition has been widely used in ID and 
2D, but in 3D several difficulties appear: not only is it computationally expensive to solve elliptic equations at each 
time step, but also one has to solve this elliptic equation with very good accuracy in order to avoid noise (see, e.g. 
flEI ). This lead to the introduction of live gauges that mimic the maximal condition near the singularity. All of the 
BM conditions are SA exactly in those cases that lead to SH formulations ||, but there are still differences between 
different subcases. For example, time harmonic slicings are not as SA [|l6| as the "I + log" are, and this is one of the 
reasons the latter has been used so much. However, if the singularity is not avoided but is excised, SA slicings are in 
principle no longer needed. On the contrary, these slicings can cause problems since they introduce steep gradients 
near the horizon. Still, there is evidence that the use of an appropriate shift can help, making the "I + log" condition 
useful even in the presence of excision |i~7| . 

Having a hyperbolic formulation with no superluminal speeds (A.; < 1) is an advantage for singularity excision since 
then one does not need to give boundary conditions at the inner boundary. In the SH families of formulations we have 
presented in this paper, all the A^'s but one (Ai) can apriori be set to 1. As we have shown, the modes associated 
with Ai are gauge modes, and thus, whether these modes are superluminal or not does not depend on the formulation 
but on the gauge condition considered. In some cases one can decide a priori whether or not Ai < I, but in principle 
Ai might depend on the solution. If that happens, one can ask the code to follow the characteristic speeds and decide 
during evolution whether boundary conditions at the inner boundary are needed or not (at the outer boundary one 
has to deal with boundary conditions in any case). If a mode does enter the computational domain and one ignores 
this fact and continues, for instance, by doing extrapolation, one is implicitly giving boundary conditions that depend 
on the grid spacing and might not have a consistent limit or be physically correct as resolution is increased. This can 
indeed happen even with very simple choices of gauge (see, e.g. the appendix of |l2]| ) and the result of such a procedure 
is uncertain. Even if one is willing to give boundary conditions to modes that enter the domain, the question of what 
conditions to give still remains. In principle, since physically relevant quantities are gauge-independent, one could 
give any boundary conditions to these modes. However, what could happen is that a bad choice leads to a gauge the 
becomes singular after a while. In fact, the same problem arises even in the absence of boundaries when the lapse is 
not prescribed apriori as a function of spacetime, but, for instance, a live gauge condition is chosen [frs) . 

An issue that we have not addressed in this paper is the introduction of dynamical shifts that may help to follow 
steep gradients near the horizon or provide some sort of minimal distortion. Current efforts are oriented along this 
line. Numerical experiments testing the formulations presented in this paper are also underway. 
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